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Abstract 

We update the study of the polarisation of Jlif/ produced in proton-proton collisions at RHIC at = 200 GeV using the QCD-based 
Colour-Singlet Model (CSM), including next-to-leading order partonic matrix elements from gluon and light quark fusion and leading- 
order conhibutions from charm-quark initiated processes. To do so, we also evaluate the corresponding cross section differential in Pj 
which agrees qualitatively with the measurements of PHENIX in the central and forward regions at low Pj ~ for instance below 2 GeV -, 
while emphasising the need for Initial State Radiation (ISR) resummation. At mid Pj, we also compare the measurements from PHENIX 
and STAR with the same evaluation complemented with the dominant or^ contributions (NNLO*). We find a reasonable agreement with 
the data. Regarding the polarisation, as shown for previous studies at larger ^fs and P^, the polarisation pattern from gluon and light quark 
fusion in the helicity frame is drastically modified at NLO and is shown to be increasingly longitudinal. The yield from charm-gluon fusion 
is found to be slightly transversally polarised. Combining both these contributions with a data-driven range for the polarisation of J I if/ from 
Xc< we eventually provide an evaluation of the polarisation of the prompt J/ifi yield which is in a good agreement with the experimental 
data from PHENIX both in the central and forward regions. 
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1. Introduction 

Until recently, the numerous puzzles in the prediction of 
quarkonium-production rates at hadron colliders were at- 
tributed to non-perturbative effects associated with channels 
in which the heavy quark and antiquark are produced in a 
colour-octet state [1, 2, 3, 4]. It is now widely accepted 
that a'l and corrections to the CSM [5] are fundamental 
for understanding the Pj spectrum of J/ ifr and T produced 
in high-energy hadron collisions [6, 7, 8, 9, 10, 11, 12]. 
The effect of QCD corrections is also manifest in the po- 
larisation predictions. While the J/ip and Y (commonly de- 
noted Q hereafter) produced inclusively or in association 
with a photon are predicted to be transversally polarised at 
LO, it has been recently emphasised that their polarisation 
at NLO is increasingly longitudinal when Pj gets larger 
[8, 10, 13, 14]. 

In a recent work [15], we have also shown that hard sub- 
processes based on colour singlet QQ configurations alone 
are sufficient to account for the observed magnitude of the 
f 7^ -integrated cross section. In particular, the predictions at 
LO [5] (Fig. 1 (a)) and NLO [6, 7, 8] (Fig. 1 (b,c,d)) ac- 
curacy are both compatible with the measurements by the 
PHENIX collaboration at RHIC [16, 17] within the present 
uncertainties. This provided some indications that the com- 
putations are carried in a proper perturbative regime. This 
agreement is improved when hard subprocesses involving 
the charm-quark distribution of the colliding protons are 



taken into consideration. These constitute part of the LO 
(al) rate (Fig. 1 (e)) and are responsible for a significant 
fraction of the observed yield, as we have argued in [15]. 

In this Letter, we proceed to the evaluation of the Pj de- 
pendence of the polarisation of the J /if/ produced at RHIC, 
both in the central and forward rapidity regions. In the sec- 
tion 2, we expose the procedure used to compute the yield 
and the polarisation at NLO (up to a*) from gg & gq fusion 
and from eg fusion at LO (at a^) as well as the procedure 
to obtain a first evaluation of the leading Pj contributions 
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Figure 1: Representative diagrams contiibuting to ■'^i 
quarkonium (denoted Q) hadroproduction in the CSM by 
gluon fusion at orders o-^ (a), a* (b,c,d) and initiated by a 
charm quark at orders (e), a* (f). The quark and anti- 
quark attached to the ellipsis are taken as on-shell and their 
relative velocity v is set to zero. 
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from gg and gq fusion at o-^ (NNLO*). In the section 3, 
we present our results. First, we show the yields differential 
in the rapidity y and Pj along with the polarisation vs Pj 
from gg & gq and from eg fusion separately. Then we show 
our results when they are combined. Finally we combine 
the latter predictions for the direct yield polarisation with 
an essentially data-driven estimation of the polarisation for 
J/ 1/' from Xc and we compare our results with the PHENIX 
data in both rapidity ranges. The last section gathers our 
conclusions. 



2. Cross-section and polarisation evaluation 

2.1. Cross section 

In the CSM [5], the matrix element to create a ^Si 
quarkonium Q of momentum P and polarisation A accom- 
panied by other partons, noted j, is the product of the am- 
plitude to create the corresponding heavy-quark pair, a spin 
projector N{A\si, S2) and R(Q), the radial wave function at 
the origin in the configuration space, obtained from the lep- 
tonic width [18], namely 



M(ab ^ Q\P) + j) ^ ^ 



N{A\suS2) 6'' RiO) 



V»^e Va^ V4^ (1) 



where P - pg + pg, p - (pg - Pg)/2, si,S2 are the heavy- 
quark spin and 6" / is the projector onto a colour- 
singlet state. In the non-relativistic limit, N{A\si, S2) can be 

written as — v(|, S2)t''m(|, ^i) where is the polari- 
sation vector of the quarkonium. The sum over the quark 
spin yields to traces evaluated in a standard way. 

In our evaluation, we use the partonic matrix elements 
from Campbell, Maltoni and Tramontano [6] to compute 
the LO and NLO cross sections from gluon-gluon and light- 
quark gluon fusion. We guide the reader to [6] for details 
concerning the derivation of the cross section at a^, the cor- 
responding expressions at o-^ can be found in [19]. In the 
case of the eg fusion (at LO), we use the framework de- 
scribed in [20] based on the tree-level matrix element gen- 
erator MADONIA [21]. 

Moreover, in order to illustrate the expected impact of 
NNLO QCD corrections for increasing Pj-, we also present 
the results when the leading-Pj- contributions at ff^ evalu- 
ated along the lines of [10] are added to the previous con- 
tributions. At this order, the last kinematically-enhanced 
topologies open up. These exhibit a P^"* fall off of their 
differential cross section in P^, typical of a single parti- 
cle exchange in the t channel. Going further in the as ex- 
pansion cannot bring any further kinematical enhancement 
as regards the Pj dependence. As a consequence, above 
a^, usual expectations for the impact of QCD corrections 
would then hold and these should be then well taken into 
account hy a K factor multiplying the yield at NNLO ac- 
curacy, which would be independent of Pj and of a similar 
magnitude as those of other QCD processes. In other words. 



a further cross-section modification between the NNLO and 
N^LO results by an order of magnitude would not be ac- 
ceptable at any Pj, while it is expected to be so between the 
cross sections at LO (Pf), NLO (P^^) and NNLO (f ^'*) for 
large enough Pj-, simply owing to their different Pj scal- 
ings. 

The procedure used here to evaluate the leading-P^ 
NNLO contributions is exactly the same as in [10]. Namely, 
the real-emission contributions at or^ are evaluated using 
MADONIA by imposing a lower bound on the invariant- 
mass of any light partons (sij). For the new channels open- 
ing up at ff^ with a leading-P^ behaviour, and which specif- 
ically interest us, the dependence on this cut is to get smaller 
for large Pj since no collinear or soft divergences can ap- 
pear there. For other channels, whose Born contribution is 
at or a^, the cut would produce logarithms of 
which are not necessarily small. Nevertheless, they can 
be factorised over their corresponding Born contribution, 
which scales as P^.^ or Pj.^, and are thence suppressed by 
at least two powers of Pj- with respect of the leading-Pj^ 
contributions (Pj'^). The sensitivity on .s™" is thus expected 
to come to nothing at large Pj. This argument has been 
checked at (NLO vs. NLO*) for T [10] and i/r [12] in the 
inclusive case as well as in association with a photon [14]. 

For the parameters entering the cross-section evaluation, 
we have taken \Rj/^(0)\^ = 1.01 GeV^ We also take 
Br(7/i/' -> f'^r) - 0.0594. Neglecting relativistic cor- 
rections, one has in the CSM, Mj/^ = 2otc. The uncer- 
tainty bands for the resulting predictions at LO and NLO 
are obtained from the combined variations of the heavy- 
quark mass within the ranges nic = 1.5 + 0.1 GeV, the 
factorisation fif and the renormalisation yU« scales chosen' 
in the couples ((0.75, 0.75); (1, 1); (1, 2); (2, 1); (2, 2)) x nir 
with nij - 4mg + P^. The band for the NNLO* is ob- 
tained for simplicity using a combined variation of m^, 
O.SniT <ijR^fif < IniT and 2.25 < < 9.00 GeVl 

2.2. Polarisation 

The polarisation parameter a (also called A) is com- 
puted by analysing the distribution of the polar angle 9 be- 
tween the {'^ direction in the quarkonium rest frame and the 
quarkonium direction in the laboratory frame. This defini- 
tion of is referred to the analysis of the polarisation (or 
spin-alignement) in the helicity frame^ . 

The normalised angular distribution /(cos 0) reads 

/(cos 0) = ^ (l+acos^0), (2) 
2(a + 3) 

from which we can extract a bin by bin in y 01 Pj. The 
interpretation of a is immediate when one relates it to the 
polarised cross sections: 

0-7- - 2o-L 



a = 



(3) 



' In principle, the renormalisation scale ambiguity can be removed us- 
ing the method described in [22]. 

^Other choices are possible and the reader is guided to [23] for more 
details. 



For ctt » cTi, thus for a yield purely transversally po- 
larised, a ^ \ while for cri » a-j (a yield purely longitu- 
dinally polarised), a - -I. On the contrary, if no direction 
is favoured, one expects ctj - {ctj^ + itt- ) = 2cr^, which 
corresponds to a ^ 0, and thus no dependence. We em- 
phasise here that the relations Eq. (2) and Eq. (3) do not 
depend on the definition chosen for 6 (the frame definition). 
Yet, the results obtained do depend on it: a yield transver- 
sally polarised in one frame can be longitudinal in another. 
In particular, the expressions that one would obtain for cr^ 
and (Tt change from one frame to another. 

3. Results 

3.1. gg and gq channels 

We first present the results from the gg and gq chan- 
nels at LO (q'^) and NLO (al+a*) in terms of differential 
cross sections as function of y and P^, which we compare 
to the experimental data for prompt J/tlr multiplied by the 
expected fraction of direct J/i// (59 + 10)% (see [15]). 

As discussed in [15], the experimental data for the Pt 
integrated cross section from PHENIX [16, 17] are quali- 
tatively well reproduced by the highest value of the theo- 
retical bands, both at LO and NLO (see Fig. 2). Let us also 
note that choosing a lower value than 1 .4 GeV for the charm 
quark mass - along the same lines as what has been done in 
studies of charm production at RHIC [25]- would lead to 
higher cross sections. However a part of the increase of 
the cross section could be attributed to the approximation 
Mj/^ = 2otc, which would then induce an artificially low 
bound-state mass in the kinematics. Such an effect is not 
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Fi gure 2: dcr jdy X Br in pp collisions at ^/sJ/n — 200 GeV 
at LO and NLO accuracy compared to the PHENIX 
data [16, 17]. The theoretical-error bands for LO and NLO 
come from combining the uncertainties resulting from the 
choice of fi f, fir, niq. We have used the LO set CTEQ6_L [24] 
for the yield at LO, the NLO set CTEQ6_M [24] for the yield 
at NLO. 



easily quantifiable without a dedicated study. Therefore, we 
do not use values lower than 1.4 GeV for here. 

Before discussing the results of the cross section differen- 
tial in Pj at LO and NLO, it is important to note that we are 
working in the low Pj region where the perturbative expan- 
sion is not always reliable. While the comparison between 
the results at LO and NLO integrated in Pj shown on Fig. 2 
gives some good indication that the perturbative expansion 
works well, it may not be so when one considers the cross 
section differential in Pj- First, the yield at NLO has shown 
negative -thence unphysical- values in some bins in Pj and 
y. When this occurs, this is always for the lowest Pj bin. 
This shows that the virtual corrections at (with a nega- 
tive relative sign) are large. This may also be the case for 
the virtual corrections at or^ which are currently unknown 
and possibly with an opposite sign to the ones at c^. This 
is a known issue which can be partly solved by resumming 
Initial State Radiation (ISR). Such a resummation has been 
carried out for the T production at the Tevatron using the 
Colour Evaporation Model [26]. One of the main outcomes 
of the latter study is that the Pj dependence of the cross 
section is significantly affected at low Pj (up to roughly 
my 12) by ISR. A similar impact of ISR is expected for the 
charmonium family, albeit for a narrower range in Pj. 

Since we have not carried out this resummation^, the 
results that we have obtained for Pj below mj/^/2 for 
dcr/dPr shown on Fig. 3 (as well as later for the polari- 
sation) are to be taken with a grain of salt. That being said, 
one sees on Fig. 3 that the slope for the yield at NLO is 
milder than for the LO, though too steep to reproduce the 
data from PHENIX [16] and STAR [27]. Yet, the general 
effect of ISR is to lower the cross section at very low Pj 
and to rise it a little up to Pj- ^ 1 - Inic. 

In addition, one expects a significant contribution from 
ff^ contributions at larger Pt- A complete evaluation of 
NNLO corrections is not yet available. For now, we can 
only rely on a study of the leading-f j- NNLO contributions 
(NNLO*) as done in [10, 12]. However, for such low val- 
ues of Pt, the latter method is a priori not reliable and we 
could not extend its evaluation (red band) in Fig. 2 (a) be- 
low Pt - 5 GeV. Overall, the comparison presented allows 
us to think that a more complete analysis (through ISR re- 
summation and matching with leading Pt contributions for 
instance) could show that the Pt dependence given by the 
CSM agrees with the experimental measurements in the low 
and mid Pt regions at RHIC energies. 

We now turn to the discussion of the polarisation results 
vs Pt (Fig. 4) in the helicity frame. As can be seen on 
Fig. 4, the complete modification of the polarisation pattern 
between the LO and the NLO results observed in previous 
works [8, 10, 13, 14] is confirmed, despite the error band 



'Such a resummation is beyond the scope of this work and will be the 
object of future investigations. 

''We did not show the NNLO* band in Fig. 3 (b) since the range in Pj 
is limited: STAR measurements are focused on the central region for the 
time being. 
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Figure 3: dcr/dPx/dy X Br from the gg and gq fusion channels in pp collisions at various order in as at ys^^ = 200 GeV: 
a) in the central (|y| < 0.35), b) in the forward (1.2 < \y\ < 2.2) regions compared to the PHENIX [16] and STAR [27] data. 
The theoretical-error bands come from combining the uncertainties resulting from the choice of fj.f, jUr, niq, see text. 
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Figure 4: aiPj) at LO and NLO from gg and gc fusion in pp at ^s^n - 200 GeV: a) in the central (\y\ < 0.35), b) in the 
forward (1.2 < \y\ < 2.2) regions. 



due to the usual mass and scale dependence on top of statis- 
tical fluctuations (see below). The LO transverse yield be- 
comes increasingly longitudinal at NLO for increasing Pj- 
In the forward region, the yield at NLO is transverse at low 
Pj and becomes longitudinal as soon as Pj > m,.. 

As discussed previously, the result obtained in the lowest 
Pj bin (below 1 GeV) are likely to be subject to higher QCD 
corrections and the band should not be extrapolated down to 
Pt - 0. Second, the extraction of the yield polarisation at 
NLO at low Pt is highly demanding in terms of computer 
time. Typically, the evaluation of the yield polarisation in 
a single Pt bin of 500 MeV with statistical fluctuation less 
than 20 % requires O(IQ^) numerical evaluations of the in- 
tegrand. These fluctuations add up to the usual theoretical 



uncertainties in the bands shown on Fig. 4. 

3.2. eg channel 

We first present the results from the eg LO (a^ ) in terms 
of differential cross sections as function of y and Pt, which 
we compare to the experimental data for prompt Jjij/ multi- 
plied by the expected fraction of direct J 1 1//. 

In [15], we have computed the yield integrated over 
Pt from eg channels and used the LO set CTEQ6.5c [28] 
based on a recent global PDF fit including Intrinsic Charm 
(IC). More precisely, we have compared the results for 
three choices for the charm distribution: (i) without IC 
[e(x,po) = (juo = 1.2 GeV)], (ii) with BHPS IC [29] 
({x}c+c= Jq x[e(x) + e(x)]dx -2%) and (iii) with sealike IC 




Figure 6: dcr/dPr/dy x Br from eg fusion in pp at -\fsMN - 200 GeV using a sealike charm distribution in the central and 
forward rapidity regions compared to the PHENIX [16] and STAR data [27] data. The theoretical-error bands come from 
combining the uncertainties resulting from the choice of fi/, yU,-, irig. 
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Figure 7: aiPj) for J /iff directly produced by eg fusion in at ^/sNN - 200 GeV using a sealike charm distribution. 



({x}c+c - 2.4%). In the following, we shall carry on our 
study of the polarisation with the sealike IC. The result are 
qualitatively the same with other choice for e(x). 

By comparing the data to our results for the contribution 
from eg fusion shown on Fig. 5, one sees that it accounts 
for 5 to 40 % of the observed yield depending on the usual 
theoretical uncertainties. Comparable fractions are obtained 
for other e{x) as discussed in [15]. 

Regarding the Pj dependence shown in Fig. 6 for both ra- 
pidity regions, one observes that the contribution is falling 
too fast compared to the PHENIX data. Nevertheless, the 
appearance of NLO contribution such as eg J/iJ^cg 
(e.g. Fig. 1 (f)) at larger^ Pj, and to a lesser extent the effect 
of ISR at small Pj, is expected to smear the curve out to 
larger Pt- 

'These were for instance analysed in the fragmentation approximation 
in [30]. 



To what concerns the polarisation pattern (Fig. 7), in the 
central region, it is similar to the one computed for J/ifr + ec 
in [7], namely nearly unpolarised as soon as Pj > nic, al- 
though not completely unpolarised with a ^ 0.2 for mid Pt- 
In the forward region, it becomes somewhat more transver- 
sally polarised and shows a dependence on the mass and 
scale choices as indicated by the widening band. 

3.3. Polarisation of the direet yield 

In order to obtain the polarisation of the direct yield, one 
has to combine the polarisation from both contributions tak- 
ing into account their proper weight. Schematically one 
has: 



(4) 



where Act is the differential cross section for one given pro- 
cess integrated in a given bin in y and Pj- 




Figure 8: dcr/dPr /dy xBr in pp collisions at ^/s^ - 200 GeV a) in the central (\y\ < 0.35) and b) in the forward (1.2 < |y| < 
2.2) regions at NLO+ (and eg LO + NNLO* for (a)) compared to the PHENIX [16] and STAR [27] data. The theoretical-error 
bands come from combining the uncertainties resulting from the choice of ///, fir, rriq, see text. 
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Figure 9: Comparison between aiPj) for J I if/ directly produced at NLO^ (see text) and the PHENIX data [31, 32] in both 
rapidity regions in pp at yfsNN - 200 GeV. 



Beforehand, we present results for the NLO^ yield, 
namely the yield at NLO accuracy from gg and gq fusion 
added to the yield from eg fusion at LO accuracy*'. The 
sum of both contributions differential in Pj is compared to 
the PHENIX and STAR data on Fig. 8 a) and b). In the cen- 
tral region, the yield at NNLO* from gg and gq fusion is 
also shown (with eg at LO added). The computation being 
close to the data for Pj < 4- GeV, it is reasonable to com- 
pare them to the PHENIX measurements [31, 32] as done 
on Fig. 9. 

Except for the lowest Pj point in the central region, 
the direct NLO^ yield polarisation is compatible with the 



*Tlie NLO corrections to eg —> J/i//X are not yet known. 



(prompt) data^, indicating a small impact of the i^(25) and 
Xc feed-downs on the polarisation. The same conclusion 
holds from the .s-channel cut analysis [33, 34], where a 
good agreement with experimental data was also obtained, 
at least in the central region. It is worth recalling that the 
latter analysis was done by neglecting the usual contribu- 
tion from the CSM which are the purpose of this work and 
are evidently not negligible. Since the results for the polari- 
sation are similar in both analysis, a combined study would 
follow the same trend. 



'Note also that for the Pt bin from 2 to 4 GeV, the yield is not perfectly 
described by the NLO^ - the deviation is slightly larger than 1 <t. Further 
contributions, such as the a^. ones, may need to be taken into account, in 
turn altering the results presented here for the polarisation. 
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Finally, the disagreement between the data and the 
present evaluation in the lowest Pj bin is not worrisome 
in view of the expected impact of higher QCD corrections 
in this bin. 



3.4. Data-driven evaluation of the feed-down effect on the 
polarisation 

For the time being, there exists no measurement of the 
direct J /iff polarisation after extraction of both the B andxc 
feed-downs. While at RHIC energy the former contributes 
less than 5 % of the full yield [35] and thence should not af- 
fect its polarisation, the Xc feed-down (up to 30-40 %) may 
impact strongly on the observed values of a. It is however 
possible to constrain its effects by using existing data from 
Hera-B on cr^,^, /o'xc^ ^'^^ instance, relying on the dominance 
of the El transition between thexc's. and the J/i//. 

Indeed, using El dominance [36], one can obtain a first 
approximation of the yield of longitudinally (transversally) 
polarised J ftp, denoted A^''"° (resp. A^fJ"') from the parent 
Xc, in terms of simple^ relations involving the polarised Xc 
yields, these are 



and 



We should emphasise that these relations neglect effects that coufd 
arise from the difference in the definition of the helicity frame of the Xc 
and that of its decay product (the J/if/ here). These effects are expected to 
vanish for large momenta for the quarkonium. Such an approximation may 
not be accurate for the smallest Pj-'s considered here. 



+ ^(r?=»Br(V2 ^'S,y) + icrl'^;'Br(V2 ^'S^y) 
+ cr\f^BreP2^'S,y). 

Let us analyse now the two otherwise extreme cases 
a™'" and a"'" : 

from Xc horn Xz 

• ff™" is reached when A^l''l=' is maximised and A^''='' 

irom Xc 

minimised. This happens when the ^P\ yield is fully 
h - Q and the ^P2 one is fully \h\ - 2. This indeed 



gives Nj - and thus a' 



from X 



= -Hi. 



a!?'" is reached when A^''"' ' is minimised and A^'' ° 

from Xz 

maximised. This happens for a ^Pi yield fully \h\ = 1 
and the one fully h - 0. This gives: 



2- P,Br(^Pi 

-H-cr3p,Br(^P2 



'Siy) 
'Siy) 
'Siy) 



(5) 



Siy) 



Recently, the Hera-B [37] and CDF [38] collaborations 
have measured the ratio of Rij = °'xci^^O(ci^J/i!'y) ^j^j^jj 

can therefore use in Eq. (5). It is however worth noticing 
that both results are somewhat different, 7? 12 = 1.0 + 0.4 
for Hera-B for < Pj- < 2 GeV at = 41.6 GeV and 
R12 = 2.5 ± 0.1 for CDF for Pj > 4.0 GeV at = 1-96 
TeV. Since we focus here on low Pj data, we prefer to opt 
for the Hera-B value, which gives for the central value 
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Figure 5: dcr/dy x Br from eg fusion in pp at -s/J^n - 
200 GeV using a sealike charm distribution compared to the 
PHENIX data [16]. The theoretical-error bands for LO and 
NLO come from combining the uncertainties resulting from 
the choice of fif, pr, niq. 
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19 



-0.47 (6) 



Supposing that about 30% of the y/t/r come from Xc de- 
cays independent of Pj in the range considered here, we ex- 
pect a partial contribution to the polarisation ranging from 
0.3 X (+1) to 0.3 X (-0.47). Regarding the other 70%, one 
can simply multiply the result obtained above for the direct 
yields (section 3.3) by 0.7, since in a good approximation 
the polarisation of Jjili from i/'(25) is expected to be identi- 
cal to the direct one. 

3.5. Comparison with the PHENIX data 

Combining the result obtained for the direct yield at 
NLO^ and the extreme expected polarisations of J lip from 
Xc, one can compare (Fig. 10) our results for the CSM 
prompt J I if/ yield with the PHENIX data in the central [31] 
and the forward [32] regions. Within the present admittedly 
large theoretical uncertainties, the agreement with the data 
is good. 
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Figure 10: Comparison between the extrapolation of a for prompt J/ifr in pp at -\/Tnn - 200 GeV (blue band), the range of 
a for the direct NLO^ (two dashed lines) and the PHENIX measurements in both rapidity regions [31, 32]. 



Three comments are in order: a) the result should not be 
straightforwardly extrapolated to Pj = 0, b) at larger Pj, 
(for gg + gq) and (for eg) will eventually dominate due 
to their Pj'^ scaling. While the former is expected from sim- 
ilar studies [10, 11, 12] to be longitudinal, we have qualita- 
tively checked that the yield from the latter (more precisely 
from eg fusion at NLO*) is similar to eg — > Jjij/e at LO; 
c) for BHPS e(x) (showing a peak of the charm distribution 
around x 0.1 - 0.3), one expects a relative enhancement 
of the yield from eg — > Jj^X for increasing Pj, i.e. for a 
momentum fraction of the charm quark in the proton ap- 
proaching the value where c^^^^(ji:) peaks. Yet, a dedicated 
study of these effects is needed to draw more quantitative 
conclusions. Indeed, depending on the relative yield of gg 
& gq fusion at NNLO* and eg fusion at LO (or NLO*) at 
large Pj, the polarisation may end up to be strongly longitu- 
dinal or slightly transversal. In any case, as opposed to non- 
relativistic QCD analyses of including colour octets [39] 
(see [40] for a recent application for RHIC energies) tran- 
sitions, the CSM yield is not expected to become strongly 
transverse for increasing Pj. 

Our results for gg and qg fusion are in qualitative 
agreement with the analysis carried out by Baranov and 
Szczurek [41] using the k, factorisation approach. However 
the latter study predicts a dominance of^c2 feed down over 
the prompt Jjij/ yield which may conflict both with the pre- 
liminary upper value [35] for the Xc feed down (< 42% at 
90% C.L.) by the PHENIX coUaboration and with the Hera- 
B results for the ratio cr^,^, l^xd- We further note that some 
contributions appearing at ff^ and enhanced by log(s) [42] 
also produce a yield which is mainly longitudinal. 

4. Conclusions and outlooks 

In conclusion, we have carried out the first NLO analysis 
in the Colour-Singlet Model of the polarisation of the direct 



Jji]/ yield at RHIC including contributions from c-quark- 
gluon fusion. Our result for the yield differential in Pj is 
in near agreement with the measurement at low and mid 
Pt both in the central and forward rapidity regions. As re- 
gards the polarisation, our evaluation for the direct yield is 
in good agreement with the PHENIX data for prompt J/t//. 
This therefore points at a small impact of the feed-downs on 
this observable in this kinematical region. 

Using constraints from existing data, namely the frac- 
tion of feed-down and the ratio cr^,^, /cr^ we have deter- 
mined likely extreme values of the polarisation of the J /iff 
from Xc by relying on the El dominance of the transition 
Xc — > J/ 4^7- This has enabled us to extrapolate our po- 
larisation evaluation of the direct yield to the prompt one. 
The results obtained are also in good agreement with the 
PHENIX data. 

Motivated by this agreement, the next step would be to 
study the Cold Nuclear Matter (CNM) effects on the yield 
polarisation in the CSM framework. These effects are ex- 
pected to come mainly from the shadowing of parton distri- 
butions and the final-state interactions between the ec pair 
and the CNM. We could for instance take benefit of the 
Glauber framework JIN [43] used in [44] to study the CNM 
effects on the J/ip yield in dA and AA collisions at RHIC; 
the latter can indeed deal with the polarisation information. 
Such a study would surely be very expediently done in order 
to refine the arguments presented in [45] to employ quarko- 
nium polarisation in heavy-ion collisions as a possible sig- 
nature of the quark-gluon plasma. 
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